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ABSTRACT 

We study stellar assembly and feedback from active galactic nuclei (AGN) around the epoch of 
peak star formation (1 <z<2), by comparing hydrodynamic simulations to rest-frame (TV-optical 
galaxy colours from the Wide Field Camera 3 (WFC3) Early-Release Science (ERS) Programme. 
Our Adaptive Mesh Refinement simulations include metal-dependent radiative cooling, star for- 
mation, kinetic outflows due to supernova explosions, and feedback from supermassive black 
holes. Our model assumes that when gas accretes onto black holes, a fraction of the energy is used 
to form either thermal winds or sub-relativistic momentum- imparting collimated jets, depending 
on the accretion rate. We find that the predicted rest-frame C/V-optical colours of galaxies in the 
model that includes AGN feedback are in broad agreement with the observed colours of the WFC3 
ERS sample at 1 <z<2. The predicted number of massive galaxies also matches well with obser- 
vations in this redshift range. However, the massive galaxies are predicted to show higher levels 
of residual star formation activity than the observational estimates, suggesting the need for further 
suppression of star formation without significantly altering the stellar mass function. We discuss 
possible improvements, involving faster stellar assembly through enhanced star formation during 
galaxy mergers while star formation at the peak epoch is still modulated by the AGN feedback. 
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1 INTRODUCTION 

The present-da y galaxy luminos ity function shows a rapid turnover at 
luminosity L* dCole et al]|200lh . while ACDM cosmol ogy predicts a 
powe r-law form for the dark matter halo mass function d Jenkins et alj 
1200 lh . The low mass-to-light ratio in low- and high-mass halos is 
interpreted as evidence for the existence of baryonic feedback, that 
regulates star formation (e.g. iBinnevI 12004). Consensus favours a 
picture in which energy from active galac tic nuclei (AGN) regu- 
lates the formation o f massive galaxies (e.g. ICiotti & Ostrikeill 19971; 
iKavirai et ai]|201 lab , despite the significant difference in the scale 
of a supermassive black hole (SMBH) from that of its host galaxy. 



* e-mail: taysun.kimm@astro.ox.ac.uk 



This picture appears to be supported by the tight observed corre- 
lation between SM BH mass and the bulge of its host galaxy (e.g. 
lHaring & Rixll2004l) . evidence for the suppression of gas cooling and 
star formatio n due to AGN interaction wit h the interstellar/intracluster 
medium (e.g. McNama ra & Nulsenl2007l) . high- velocity galactic out- 
flows that are commonly de tected in quasars dPounds et al.l [2003; 
iGangulv & Brothertonl l2008h and X-ray cavities inflated by radio 
jets, associated w ith AGN, that are observed in cluster centres 
jFabian et alj2006t) . Nevertheless, it is still unclear as to how, to what 
extent, and at what epochs, AGN feedback shapes the star formation 
history (SFH) of galaxies in the Universe. 

Several studies have used semi-analytic or fully hydrodynamic 
models of galaxy formation, to u nderstand the effect of negative 
feedback on the cosmic SFH (e.g. iFontanot et alj|2009l ). While de- 
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tailed physical ingredients and prescriptions differ in each model, 
such studies agree on the importance of AGN feedback for repro- 
ducing the masses, co lours, and star form ation rates (SFRs) of local 
massive galaxies (e.g. lKimm et alj[2009h . However, the comparisons 
are almost exclusively performed in the nearby Universe, while the 
bulk of the stellar assembly in these galaxies takes place at high red- 
shift (z > 1). Unfortunately, 8-10 Gyrs of evolution can easily wash 
out many of the details of stellar assembly, and thus it is neces- 
sary to investigate the epoch at which the bulk of the star formation 
takes place (I <z<3) dMadau et al Jl 19981) . Confronting the observed 
galaxy colours with the models at these epochs represents a more 
stringent test of their reliability. 

In this Letter, we study the stellar assembly in massive galax- 
ies, by comparing a fully hydrodynamical cosmological simulation 
to rest-frame UV-optical galaxy c olours from the WFC3 Early Re- 
lease Science (ERS) programme (iWindhorst et al.lfeoi ll) . The rest- 
frame UV (shortward of 3000A), which is sensitive to even resid- 
ual amounts of star formation, is a powerful tracer of how q uies- 
cent a galaxy is (e.g. lKavirai et all2011bl ; lRutkowski et alj2012t) . en- 
abling us to constrain the quenching of star formation due to AGN 
and stellar feedback. By contrast, the rest-frame optical colours con- 
strain the average epoch of stellar mass assembly. This is the first 
direct comparison to simulations of rest-frame UV-optical colours 
at the epoch of peak star formation. We describe the data and sim- 
ulations in § [2] and [3] respectively. We present the comparison be- 
tween observations and models in § [4] and discuss the implica- 
tions of our results on galaxy formation in § [5] Throughout, we 
adopt (ho, a,„, £l b , Q a, or) = (0.7, 0-26, -045, 0.74, 0.8) following 



the WMAP-3 results fainshaw etai]|2009i) . All fluxes are based on 
the AB magnitude system. 



2 OBSERVATIONS 



The WFC3 ERS programme dWindhorst et al.l201 lh has imaged ~ 1/3 
of the GOODS-South field (~45 arcmin 2 ), using the WFC3 F225W, 
F275W, F336W, F098M, F125W, and F160W filters, with expo- 
sure times of 1-2 orbits per filter. In combination with the exist- 
ing GOODS BViz coverage dGiavalisco et alj [20041) . this provides 
10-filter panchromatic coverage (0.2-1.7 pm) to point source depths 
of AB -26.5-27.5 mag (UV-IR, 5a). Here, we study 188 and 151 
galaxies with M star > W w Mr. } , for comparisons at 0.8 < z < 1.2 and 
1.7 < z < 2.3 respectively. Note that the ERS sample is complete 
down to ~ 10 9 Mq and the observed images that trace both the rest 
frame UV and the rest frame optical are deep enough for all massive 
galaxies to be detected (Windhorst et al.l201 ll) . For every galaxy, pho- 
tometric redshifts were calculated using EAZY, based on all 10 filters 
with all possible com binations of pairs of the EAZYvl.O templates 
terammer et al.ll2008h . The absorption due to intervening intergalac- 
tic HI clouds is taken into consideration. No redshift prior probability 
distribution specific to these fields is used. The resulting mean re- 
duced x 2 °f the fit is 0.94, and the typical redshift uncertainty of the 
sample is Az = |z spe c — Zphot | — 0. 15. We note that no k-correction is 
applied to the observed sample, and that we refer the rest-frame UV 
and optical fluxes to the fluxes through the WFC3/ACS bandpass that 
roughly trace the GALEX NUV or Johnson B/V. Thus our estimated 
rest-frame NUV or optical fluxes are slightly different from GALEX 
NUV or Johnson B/V fluxes. 

In addition, the full suite of WFC3/ACS photometry of each 
galaxy is compared to a library of exponentially decaying m odel 
star formation histories, based on the lBruzual & Chariot! d2003h stel- 
lar models. A wide range of ages (0.001-13 Gyr), decay timescales 
(0.1-9 Gyr), metallicities (0.005-2 Z s ) and dust extinctions < 



E(B — V) < 2 are employed, with the normalisation of the models 
yielding the stellar mass. The likelihood of individual models (% 2 /2) 
are calculated and parameters are marginalised. The median of the 
marginalised probability distributions are taken to be the best esti- 
mates, with the 25 and 75 percentile values (which enclose 50% of 
the total probability ) providing an associated uncertainty (see e.g. 
iKavirai et alj|201 lbl) . Specific SFRs (SSFRs) and stellar masses thus 
derived have uncertainties of ~0.3 dex. 



3 SIMULATION 

In this paper, we use the oc tree-based Eulerian hydrodynamics 
code, Ramses dTevssier 2002). Metal-dependent radiative cooling 



is modelled based on 



Sutherland & Dopital d 19931) . and a uniform 
UV background is i n stanta neously turned on at z = 10.5 follow- 
ing |HaardT&Midau| l ll996l) . Star particles are created following a 
Schmidt law with 2% efficiency, when the density of a gas cell 
reaches a critical density, po = 0.4 H/cm 3 . We use the polytropic 
equation of state for the gas above the threshold density, with the min- 
imum temperature of 10 4 K. Massive stars are assumed to lose their 
mass through stellar winds and a supernova phase, dispersing gas 
and metals (18% and 2% by mass, respectively) into their surround- 
ings. Feedback is modelled as an isotropic kinetic outflow, which 
carries matter amounting to ten times the mass loss from the stars 
toubois & Tevssiej|2008l) . 

Our model assumes that seed BHs of 10 5 M S form in re- 
gions o f high gas density, ensuring that only one BH is formed per 
galaxy jDubois et al.1 120121) . The growth of the BH is tracked self- 
consistently, based on a modified Bondi-Holye-Lyttleton accretion 
rate. We adopt a density - depen dent boost factor ocbh = (p/Po) 2 , fol- 
lowing HoothAlch^i j2009h . When gas accretes onto BHs, we as- 
sume that a central BH can impact the ambient gas in two ways, de- 
pending on the Eddington ratio (% = m/ih e ^). For a high accretion 
rate (% > 0.01), 1.5% of the accretion energy is injected as thermal 
energy (quasar-like mode), whereas sub-relativistic (v = 10 4 km/s) 
momentum-imparting collimated winds are laun ched for a low accre - 
tion rate (% < 0.01) event with 10% efficiency dDubois et aTll2012h . 
The parameters are chosen to match the stellar baryon fraction in 
groups and clusters of galaxies (Dubois et al. in prep). In order to 
avoid artificial cooling and maximise the impact of AGN feedback, 
we store the thermal energy and kinetic energy until it either reaches 
10 8 K or the energy corresponding to 10% of the BH mass. There- 
fore, the efficiency of the two modes should be regarded as a rough 
estimate, not a true efficiency. 

We have performed two cosmological simulations with th e same 
initial conditions generated by mpqr afic jPr unet et al. 2008), a par- 
allel version of the graf ic package dBertschingej200 lh . The simula- 
tions contain 256 3 dark matter particles with m^ m ~ 4 x 10 9 M© in a 
50 hT x Mpc comoving periodic box. We fix the maximum spatial res- 
olution to ~ 0.38 kpc/h in physical units. The minimum mass of the 
star particles is 1.65 x 10 6 M R . The virial mass of the most massive 
halo in our simulated volume is 3.3 x lO'^M© at z = 1. Note that the 
mass of dark matter halos used to compare with the observation in §4 
is in the range 12< logM v ; r < 13.5, and hence our sample represents 
massive galaxies (m stai > 10 M Q atz = 1) forming in group/field en- 
vironments. This allows for a fair compariso n with the ERS sample , 
which is known to have few cluster galaxies dWindhorst et al 
It is also worth mentioning that 36% of the simulated massive galax- 
ies are bulge-dominated at z = 1 , which seems compat ible with the 
obser ved morphological fraction in the GOODS fields dBundv et al.1 
l2005h . 
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( F850LP - F098M ) ( F125W - F160W ) 

Figure 1. Colour comparisons of massive galaxies at z = 1 (left; m stal - > 10 M@) and z = 2 (right; m siM > 10 5 M ) with simulations that include AGN+supemova 
feedback. The observed colours on the x and y axes trace rest-frame (B — V) and (NUV — V), respectively (see text for more details). Filled circles with error bars 
indicate the WFC3 ERS data, while simulated galaxies are shown using filled grey circles. The size of the symbol scales with stellar mass. The three lines represent 
the expected colours of simple stellar populations (SSPs) with different ages at different redshift: in the left-hand panel, z =0.8 (dashed), 1 .0 (solid), and 1 .2 (dotted), 
in the right-hand panel: z = 1.7 (dashed), 2.0 (solid), and 2.3 (dotted). For clarity, the colours of SSPs with different ages (0.1, 0.2, 1, and 4 Gyr) are indicated by 
filled square symbols. Points with a concentric circle denote galaxies with low specific star formation rate (SSFR< 0.01 Gyr -1 ). We indicate the median uncertainty 
of the observed colours in the top right and the typical extinction vectors of the simulated galaxies in the left of each panel. 



Predicted colours of simulated galaxies are generated by 
folding the stellar spectra o f individual star particles with the 
stellar population models of iBruzual & Charloa d2003l) . The op- 
tic al depth is calculated following the empirical calibration 



of iGuiderdoni & Rocca-Volmerangel 1 1987[). an d dust attenuation 



calculated using the ICardelli et al.l 41989b extinction 



jPevriendt et al.l2O10L Eq. 1). We emphasise that the attenuated spec- 
tra of the simulated galaxies are redshifted to either z = 1 or z = 2 
before convolving them with the throughput of WFC3/ACS filters to 
directly compare with ERS galaxies. 

4 RESULTS 

In Fig. [TJ we present the rest-frame UV and optical colours of ERS 
galaxies around z ~ 1. We use the ACS F435W, F850LP, and WFC3 
F098M filters to trace the rest-frame NUV (~2300A), B, and V bands, 
respectively. In order to test the influence of AGN feedback, we re- 
strict our sample to galaxies more massive than m s tar > 10 Mg be- 
cause SN feedback is ineffective in such large potential wells, leaving 
AGN feedback as the most plausible mechanism for modulating star 
formation. Three reference lines are included in the figure to help 
distinguish the age and redshift sequence. Each line corresponds to 
a simple stellar population (SSP) with various ages (marked next to 
the line). Differences between the lines indicate the change in the 
colours when a spectrum is redshifted to different values before it 
is convolved with the filter throughputs (z = 0.8, 1.0, and 1.2 [left 
panel] or z = 1.7, 2.0, and 2.3 [right panel]). Observed galaxies show 
a large scatter in (NUV — V) and (B — V) colours, partly because no 
k-correction is applied. When the redshift interval of the sample is 
constrained around z — 1, a colour sequence is apparent (red points), 
with massive galaxies being redder in both UV and optical colours. 

To make a detailed comparison, we study the SSFRs derived 
from SED fitting (see Section 2), marking 'quiescent' (SSFR = 
'«star/"2star < 0.01 Gyr -1 ) galaxies in our sample with a concen- 
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Figure 2. The fraction of passive galaxies for the observations (black squares), 
the model with supernova and AGN (red stars), and the model with super- 
nova (blue diamonds). Top panels show the fraction of galaxies showing 
SSFR< 0.1 Gyr -1 , while bottom panels display the fraction of galaxies which 
are an order of magnitude less active SSFR< 0.01 Gyr -1 . Poisson errors are 
indicated by error bars. 



trie circle in Fig. [TJ We find that almost all massive galaxies with 
(NUV — V) >4 mag are quiescent. It is worth emphasising here that 
this provides a stringent test for any cosmological galaxy formation 
model, since such low S SFRs are comparab le to those of nearby 'red 
and dead' galaxies (e.g. iKimm et alj|2009l) . Thus, any process that 
shapes the star formation history of massive galaxies must be power- 
ful enough to quench star formation at this early epoch. 

We find that the simulated galaxies with AGN feedback agree 
reasonably well with the observed sequence at z = 1 . The left-hand 
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panel of Fig. Q] indicates that a majority of model galaxies with 
(NUV — V) > 4 mag are indeed quiescent (double circles), but not due 
to dust reddening, indicating that AGN feedback is a plausible expla- 
nation for the emergence of quiescent galaxies at z — 1 . However, in 
the absence of AGN feedback almost all model galaxies are predicted 
to be actively star forming (Fig. [2}, suggesting that gravita tional heat- 
ing by substructure alone (e.g. lKhochfar & Ostrike rll2008l) is not able 
to account for the build-up of qui escent galaxies in grou p/field en- 
vironments at this epoch (see also ljohansson et aLll2012l) . Note that 
such interactions between infall of substructure and the interstel- 
lar/intracluster medium is automatically accounted for in our hydro- 
dynamic calculations. We also see that the model with AGN feed- 
back reproduces the observed trend of blue galaxies ( (NUV — V) < 3 
mag) being less massive. However, there is an indication that massive 
galaxies with intermediate colours ((NUV — V) ~ 3.5 mag) are more 
frequent than observed in the data (see the left panel in Fig. [I). 

In Fig. Q] (right panel), we show the corresponding colour com- 
parisons at z — 2. We only show galaxies with masses greater than 
3 x 10 M@, which form more than 90% of the progenitors of our 
model galaxy sample at z = 1 . Since the redshift of the galaxy sam- 
ple has changed, we now use the ACS F606W and the WFC3 F125W 
and F160W filters to trace the rest-frame NUV, B, and V bands re- 
spectively. We find that our simulated galaxies at z — 2 compare bet- 
ter with WFC3 galaxies at 1.7 < z < 1.9 rather than 1.9 < z < 2.1. 
We suspect that this is due, in part, to a systematic u nderestimation 
of the photometric redshifts at these epochs (see also IWilliams et alj 
120091) . Note that the observed galaxies at 1.7 < z < 1.9 are compa- 
rable to the model galaxies at z = 2, given the typical uncertainty 
of EAZY redshifts at these epochs (Az ~ 0.15). At z — 2 we find 
that there exists a non-negligible number of massive galaxies with 
a very low level of star formation activity, which are absent in the 
run with AGN feedback. The massive galaxies (m sta r^ 10 11 Mq) in 
the model are actively star forming (SSFR~ 0.3 Gyr ), leading to 
colours ((NUV — V) < 3) that are inconsistent with the observations. 

Fig.|2]shows the fraction of passive galaxies, which confirms that 
the modelled recent star formation needs to be further suppressed in 
massive galaxies at z — 2. While our simulation predicts that < 20 % 
of massive galaxies have SSFR< 0. 1 Gyr -1 at z — 2, the observations 
suggest that 30-60% of the massive galaxy population is passive (top 
right). At z — 1, a majority (60-80%) of the massive galaxies in the 
model have become passive (top left), but when a more conservative 
definition for a passive galaxy is used (SSFR< 0.01 Gyr -1 ), we find 
that the passive fraction in the model drops to approximately 0-20% 
(bottom left). This is inconsistent with the observations, which sug- 
gest that 50-100% of massive galaxies are quiescent at this epoch. 

Nevertheless, the impact of AGN feedback becomes clearly vis- 
ible in the stellar mass functions (see Fig. |3}. The number of mas- 
sive galaxies in the AGN run (m sta i->10 M©; solid lines) are sys- 
tematically smaller than the run without AGN (dashed lines) by a 
factor of ~ 3. Comparing twin halos that share position and mass 
in the feedback/no feedback runs, we find that the stellar masses 
of the galaxies they host are reduced by a factor of 2-5, meaning 
that the mass functions are shifted horizontally (to the left) due to 
AGN feedback. The mass functions are in reasonable agreement with 
the stellar mass function estimated by IPerez-Gonzalez etal, I d2008h 
(red and blue shadings), and also compatible with the simple num- 
ber counts from the WFC3 galaxies at 1 .7 < z < 2.3 (orange stars) or 
0.8 <z< 1.2 (blue stars). 

One may wonder whether our AGN feedback recipe is too 
strong, given the apparent deficit of the most massive galaxies. How- 
ever, this can be partly attributed to the absence of density fluctua- 
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Figure 3. The stellar mass functions of galaxies at z=2 (reddish colours) and 
z=l (bluish colours). Our estimations at each redshift intervals for simulated 
galaxies with and without AGN feedback are plotted as solid and dashed lines, 
respectively. Error bars indicate Poisson errors in the counts. For comparison, 
we display the mas s functions at 0.8 < z < 10 (b lue shading) and 1 .6 < z < 2.0 
(red shading) fromlPerez-Gonz alez et al li liooi . Also included as stars are the 
simple number counts from the WFC3 sample. 
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Figure 4. The median star formation history of simulated galaxies at z = 1 
in the ran with (solid) and without (dotted) AGN feedback. Top: the median 
star formation rates of galaxies in three different mass bins. For clarity, the 
SFRs are scaled down by 100 for galaxies with 10" 5 < m star < 10 12 2 (red), 
50 for galaxies with 10 110 < m stiir < 10" 5 (cyan), and 20 for galaxies with 
10 10 < m star < 10" (blue). Bottom: The specific star formation rates (SSFRs) 
as a function of look-back time. 

tions on scales com parable to the size of the simulation (jsirko 2005 ; 
Gned in et al Moreover, the fact that we need further quench- 

ing of residual star formation in massive galaxies does not appear to 
favour moderating the impact of AGN feedback. We discuss this issue 

further in the next section. 

Finally, in agreement with the downsizing picture dCowie et al.l 

Il996h . we find that in the run with AGN feedback, the star formation 
histories of more massive galaxies are peaked at higher redshift (solid 
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lines in the upper panel of Fig. [4). This can also be seen in the bottom 
panel, where we plot the median SSFR of 1 127, 117, and 29 galaxies 
in three mass bins, 10 10 <m stal -< 10 n M s , 10 11 <m 8tar < 1O IL5 M , 
and 10 115 < m star < 10 12 ' 2 Mq, as a function of look-back time. If we 
compare the galaxy SFHs in the runs with (solid lines in the top panel) 
and without AGN feedback (dotted lines), we see that star formation 
is preferentially suppr essed at 1 <z<3 a t which the bulk of the star 
formation takes place dMadau et al.ll 1998h . 

Interestingly, we see that the most massive galaxies have an ex- 
tended SFH, which may, at first sight, seem inc onsistent with the hig h 
observed [oc/Fe] ratios in local systems (e.g. iThomas et al.l 120051) . 
However, we note that it is the star formation timescale in the sub- 
units that finally merge to form the larger galaxy that determines its 
average [oc/Fe]. In other words it is not possible to directly infer the 
[oc/Fe] from a 'composite' SFH of the various subunits (such as the 
one shown in Fig. @, because much of the broadness of the SFH 
is driven by the fact that star formation peaks at different epochs in 
different subunits. Hence, the degree of the discrepancy may not be 
substantial as it appears, and will be tested with more detailed chemo- 
hydrodynamic simulations in a forthcoming paper. 

5 DISCUSSION AND CONCLUSIONS 

We have presented the first direct comparison of the rest-frame UV- 
optical colours of massive galaxies at l^z^2 to the predictions 
of hydrodynamic simulations that include AGN feedback. We have 
shown that massive galaxies in the simulations are actively star form- 
ing at this redshift in the absence of a strong source of energetic 
feedback. In an attempt to estimate the stellar assembly of massive 
galax ies, we have incorpo rated a simple, jet/quasar-mode AGN feed- 
back dPubois et ai]|2012h . and found that it reproduces a reasonable 
fraction of moderately quiescent galaxies in both rest-frame UV and 
optical colours at z = 1 . However, the model slightly overpredicts the 
level of recent star formation activity in massive galaxies at z = 1 and 
2, compared with the observations. Despite this, the stellar mass func- 
tions in the model with AGN feedback agree reasonably well with 
the WFC3 ERS data (within the statistical uncertainties), due to the 
efficient quenching of star formation during the epoch of peak star 
formation. 

The lack of massive, quiescent galaxies in the model seems to 
require further suppression of residual star formation without sig- 
nificantly altering the stellar mass function. This suggests that the 
peak epoch of the stellar assembly in massive galaxies should com- 
mence earlier than we have predicted, while the residual star forma- 
tion must decline more sharply over time. Although further star for- 
mation quenching is attainable by introducing more enhanced feed- 
back or other processes that we have ignored (e.g. photoionisation of 
gas coolants due to soft X-ray and extreme UV photons from stars, 
ICantaluDoll20ld) . changing the peak epoch of stellar assembly, while 
keeping the stellar mass function fixed may require well-balanced in- 
terplay between feedback and star formation. 

Regarding the extended nature of the predicted SFHs, a differ- 
ent approach, involving a h igher level of sta r formation activity in 
the very early universe (e.g. lStarketaLll2009l) is also appealing. The 
advantages of this idea are 1) it does not violate the stellar mass func- 
tions at z < 2, 2) it alleviates the discrepancy in [oc/F e] ratio in mas- 
sive galaxies to some degree (e.g. iPipino et al .1120091) . although other 
routes, su ch as the use of a top- heavy IMF may offer more promising 
solutions ( lArrigonietai1l2010l) ). and 3) star formation can be effi- 
ciently suppressed with the same amount of energy released from the 
black holes by consuming gas earlier on. Even though the physical 
origin of such high specific star formation rates is unclear, we note 



that proper modelling of gas fragmentation in galax y mergers can 
easily elevate star formation by an order of magnitude dTevssier et al.l 
l201Ch . and may eve n be able to match the o bserved high SSFRs in the 
early universe (e.g. lKhochfar&Silkll201ll) . 

We also note that the mass and spatial resolution of a simulation 
plays a role in determining the star formation history of the model 
galaxies. Adopting a better resolution allows us to resolve smaller ha- 
los in the early universe (z>5), leading to earlier star formation, with 
an order of magnitude enhancemen t in the star formation activity at 
this epoch dRasera & Tevssiedl2006t) . This will naturally increase the 
mean age of the massive galaxies (e.g. I Johansson et al]|2012h . How- 
ever, the question is how many stars form during this epoch, compared 
to the total stellar mass at z — 1 or 2. Inspection of the top panel of 
Fig. E] shows that an order of magnitude increase in the SFR at z > 5 
will not significantly affect the overall SFHs. 

Given the simplicity of our modelling of AGN feedback, the 
agreement between observations and models is encouraging. An im- 
portant challenge will be to match chemical abundance ratios in mas- 
sive galaxies, while producing the correct number of red and dead 
galaxies at high redshift. In a forthcoming paper we will investigate 
these issues in detail, focussing on more realistic modelling of the 
starburst phases and the interplay between feedback and star forma- 
tion at high redshift. 
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